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ABSTRACT 



Aims. We probe the stellar population age and metallicity distributions in nearby elliptical galaxies over the largest extension to date. 
Methods. Long-slit spectroscopy is made by using the spectrograph SCORPIO of the 6-m telescope of the Special Astrophysical 
Observatory of the Russian Academy of Sciencies. The Lick indices Hy6, Mg b, Fe 5270, and Fe 5335 are calculated along the slit up 
to radii of 1.3 to 3 in 4 galaxies and up to 0.5 in the fifth one. The comparison with evolutionary synthesis models of simple 
stellar populations allows us to disentangle age and metallicity and to measure both. 

Results. We have found that the mean stellar age is constant along the radius only in one galaxy out of 5. The other 4 galaxies 
demonstrate quite different behaviour of the mean stellar age: the outer parts are older than the centres in 3 cases and younger - in 
one case. The metallicity gradients cannot be approximated by a single power law over the full radial extension in 4 galaxies of 5. The 
inner metallicity gradients within 0.5 are all rather steep, steeper than -0.4 metallicity dex per radius dex, and are inconsistent with 
the origin of the elliptical galaxies by a major merger. 

Keywords. Galaxies: individual: NGC3348, NGC2634, NGC2810, NGC4339, IC 2179 - Galaxies: evolution 



1. Introduction 



Elliptical galaxies were described by iHubbld (11926') as smooth 
homogeneous stellar systems. However even early photomet- 
ric studies revealed the exi stence of colour gradients along the 
radius in bright ellipticals ( 'Tifftlll969l: iBurkhead & Kalinowskil 
1 1974), which had been treated later as metallicity gradients be- 
cause of consistency between the meta llicity slopes determined 
from different-band colour grad ients jSandage & VisvanathanI 
119781: iPeletier & ValentiinI Il989h . But there were also some 
hints, after comparing colour gradients to absorption-line 
strength gradients, that there may be also age variations along 
the radius that make a minor contribution to the c olour gradients 
jGorgas. Efstathiou. & Aragon Salamancalll990h . 

The presence of metallicity gradients along the radius in el- 
liptical galaxies is a crucial test for all models of galaxy forma- 
tion. The model of mono lithic dissipati ve collapse of a proto- 
galactic gaseous cloud by iLarsoiil d 19741) predicted the metallic- 
ity gradients with the correct sign: during the violent star for- 
mation phase over the whole galaxy, the gas enriched by fresh 
heavy elements fell into the centre and became the material for 
subsequent star fo rmation. The mo re elaborated models of dissi- 
pative collapse by ICarlber3 (1 1 984t) gave the characteristic value 
for the final stellar metallicity gradient in the largest elliptical 
galaxies of -0.5 dex per radius dex, and the dependence of the 
gradient slope on the galactic mass. Sm all elliptical gala xies had 
no metallicity gradient in the models of lCarlberS ( |1984|) . due to 
the early occurrence of a galactic wind that stopped the galaxy 
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formation. The other main scenario in which elliptical galaxies 
are formed is the framework of hierarchical merging. In simula- 
tions of elliptical galaxies forming by major merging, metallic- 
ity gradients are also predicted. In general, the gradients in these 
simulations are sh allower than t hose formed i n the monolithic 
collap se scenario (IWhitdll980l: lBarneslll996l: iBekki & Shiovi 
\199% . 

Recently iKobayashil ( |2004|) has considered the evolution 
of metallicity gradients in the framework of a modern LCDM 
model of elliptical galaxy formation. The processes considered 
include monolithic collapse at high redshift and subsequent mi- 
nor and major mergers followed in some cases by induced star- 
formation bursts. Monte Carlo simulations have given 124 var- 
ious histories of elliptical galaxy evolution and reproduced the 
observed variety of metallicity gradients. The extreme case of a 
galaxy formed by the pure monolithic collapse and then evolving 
quietly corresponds to a very steep metallicity gradient, about 
A[Fe/H]/Alogr -1. During the evolution, the metallicity gra- 
dients become shallower when galaxies merge; no galaxy hav- 
ing experienced some major merger has the metallicity gradient 
steeper than -0.35 d ex per radius dex . Among the galaxy en- 
semble simulated by IKobayashil (12004 ) the average metallicity 
gradient is A[Fe/H]/Alogr ^ -0.3, with a scatter of +0.2 dex 
- quite consistent with the observations of colour gradients in 
elliptical galaxies; no correlation between gradient and galaxy 
mass is found. These simulations demonstrate that the knowl- 
edge of exact values of metallicity gradients in ellipti cal galaxies 
i s imp ortant for probing their evolution. In particular. IKobayashil 
(l2004l) suggests that the metallicity gradient in any individual 
elliptical galaxy is an indicator of its history and may reveal if 
major mergers were significant during its evolution. 

The current observational data provide mean metal- 
licity gradients A[Fe/H]/Alogr typically in the range 
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-0.2 to -0.3 (see e.g . rCarollo. Danziger. & BusonI 
Davies. Sadler. & Peletieil 119931; 'Kobava shi & Arimotol 
Mehlert et al. 2003), while observational evidence of the 
predicted correlation between metallicity gradients and galaxy 
mass or l uminosity remains somewhat controversial ( jVader et al.l 
1988; P eletier et al.lll990l;lFbrbes. Sanchez-Blazquez . & Proctor 
'2005ft . 



19931: Table 2. Long-slit spectroscopy of the galaxies 



Unfortunately, investigating the metallicity gradients in el- 
liptical galaxies is not straightforward, and these observa- 
tional results should be interpreted with caution. Indeed, most 
estimates are made under the a priori assumption that the 
whole colour or absorption-line gradient results from metal- 
licity gradient. For example , we could refer to the paper by 
iKobavashi & Arimotol d 19991) who summarise a totality of data 
on metal-line index gradients in 80 elliptical galaxies over 20 
independent studies and calculate metallicity gradients by using 
evolutionary synthesis of simple stellar populations (SSP) and 
by assuming that the stellar populations at any radius in all 80 
galaxies under consideration are 17 Gyr old and no age gradient 
exists. But we know that age-metallicity degeneracy concerns all 
the optical colours, as well as all the metal-line indices. If one 
diminishes either metallicity under a constant age or age under 
a constant metallicity, all the colours would be bluer and all the 
metal lines would be weaker. So ages and metallicity must be de- 
termined simultaneously by confronting for examples metal-line 
indices to the Balmer lines or to the broad-band optical colours, 
to make definite conclusions about the metallicity (and age) gra- 
dients. Secondly, due to the different scales one must be careful 
when comparing observationally-determined metallicity gradi- 
ents to the results of simulation. For example, the metallicity 
gradients presented by Kobavashi (2004) cover the region out 
to 2 ^e, whereas most observational metallicity gradients are de- 
rived over a relatively compact region, typically within 0.5 to 
1 re. Simple extrapolations of the measured gradients to larger 
radii are very dangerous possibly leading to erroneous conclu- 
sions. 

This work presents a precise determination of metallicity and 
age gradients in five elliptical galaxies made by using Lick- index 
measurements and by breaking the age-metallicity degeneracy 
with the two-index diagrams based on the evolutionary synthe- 
sis of old SSPs. The other feature of this work is a significant 
extension of the index profiles: for 4 of 5 galaxies we reach 1.3 
to 3 effective radii, whereas most works up to now have consid- 
ered index gradients within one effective radius. In Sect. 2 we 
present our sample and in Sect. 3 we discuss the observations, 
data reduction, and the analysis of the spectra. In Sect. 4 we dis- 
cuss the raw line-index profiles for the galaxies in our sample 
and convert these indices to age and metallicity gradients using 
SSP diagrams. These results are discussed in Sect. 5. Section 6 
offers a summary. 

2. Sample 

This work is part of the larger project involving complex pho- 
tometric and spectral studies of nearby round elliptical galaxies. 
The ultimate aim is to probe the dynamical and stellar structures 
of elliptical galaxies up to a few effective radii from the centre. 
Round ellipticals were selected because of their probably spher- 
ical shape so the simplest one. For this reason we suggest that 
one longslit cross-section with an arbitrary orientation of the slit 
would be enough to obtain a complete picture of radial varia- 
tions for all kinematical and stellar population properties. The 
database HYPERLEDA has been used to select all EO galax- 
ies with a visible ellipticity less than 0.075 (at the 25th blue 



Date 


Galaxy 


Exposure 


PA (slit) 


Seeing 






(min) 


(deg) 


(arcsec) 


6 Mar 05 


NGC 3348 


100 


34 


3 


7 Mar 05 


IC2179 


220 


38 


3 


7 Mar 05 


NGC 2810 


120 


63 


3 


20 Apr 06 


NGC 3348 


120 


293 


2 


21 Apr 06 


NGC 4339 


40 


241 


2 


26 Apr 06 


NGC 2634 


280 


45 


4 



isophote) and with the diameters, also at the 25th blue isophote, 
between 1 and 3 arcminutes. The latter condition is imposed to 
provide sufficient spatial resolution, on one hand, and to match 
the field of view of the spectrograph SCORPIO, on the other . We 
must have blank sky measurements for the external parts of the 
slit to subtract the sky contribution properly. The northern part 
of the sample lists 18 galaxies. Now we present the first part of 
the results of the spectral observations for 5 of them. The global 
characteristics of the considered galaxies are given in Table 1 . 

3. Observations and data reduction 

Long-slit spectral observations were made with the multi- 
mode focal reducer SCORPIO (lAfanasiev & MoiseevI 
200?) installed at the prime focus of the BTA 6-m tele- 
scope at the Special Astrophysical Observatory. A de- 
scription of the SCORPIO instrument can be found at 
http://www.sao.ru/hq/moisav/scorpio/scorpio.html. For our 
observations with the SCORPIO we used the phase-volume 
holographic grating (the grism) VPHG2300 providing the 
narrow spectral range of 4800-5540 A; this wavelength region 
is rich in stellar absorption lines and contains the Lick indices 
HyS, Mgb, Fe5270, and Fe5335 suitable for studying stellar 
populations. Moreover, the emission lines [Olll] AA4959, 5007 
appeared to be strong in some cases; they are used for the gas 
kinematics study and for estimating the emission contamination 
level for the Lick index H/3. The slit width was one arcsec, 
and the spectral resolution about 2.2 A. The CCD 2k x 2k 
was used as a detector, and the scale across the slit was 0.35" 
per pixel. The full log of observations is given in the Table 2. 
B esides t he galaxies, several G- and K-giant stars from the list 
of lWorth ev (1994) were observed in the same mode, and their 
spectra were used for the Lick system calibration and for cross- 
correlation with the spectra of the galaxies. The preliminary 
data reduction was performe d by means of IDL-based software 
(lAfanasiev & Moiseevll2005l) . 

During the further reduction we applied logarithmic bin- 
ning along the slit to provide a sufficient signal-to-noise ratio 
(S/N > 25-30 per bin). For every galaxy we calculated radial de- 
pendencies of the Lick indices H/3, Mgb, Fe5270, and Fe5335 
that are suitable to determine metallicity, age, and Mg/Fe ra- 
tio of old stellar populations (Worthey et al. 1994). To calibrate 
the SCORPIO index system o nto the standard Lick one, we ob- 
served 22 stars from the list of lWorthev et al.l d 19941) during sev- 
eral observational runs with the same grating and spectrograph 
setup; the comparison of the instrumental and standard Lick in- 
dices for these stars from the work of Worthev et al. (1994) is 
shown in Fig. 1. We calculated the linear regression formulae 
to transform our index measurements into the Lick system. The 
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unit 


source 


NGC 2634 


NGC 2810 


NGC 3348 


NGC 4339 


1C2179 


Type 




NED' 


El: 


E 


EO 


EO 


El 


Vr 


km s"' 


NED 


2258 


3572 


2837 


1289 


4336 


Distance 


Mpc 


LEDA^ 


34.7 


53.5 


44.3 


18.2( Virgo) 


64.6 




arcmin 


LEDA 


1.7 


1.7 


2.0 


24 


1.1 




kpc 


LEDA 


17.2 


26.5 


25.8 


12.7 


20.7 




mag 


RC33 


12.71 


13.07 


11.84 


12.31 


13.21 


Mb 


mag 


LEDA 


-20.0 


-20.7 


-21.7 


-19.3 


-20.9 


(B - V)« 


mag 


RC3 


0.90 


0.95 


0.95 


0.90 


0.93 


W - s)« 


mag 


RC3 


0.48 




0.44 


0.50 


0.52 


0-0 


km s"' 


LEDA 


182 + 4 


223 ±8 


238 ± 10 


113 ±4 


146 ± 52 




arc sec 


(see notes) 


18'* 


145 


22-* 


26' 


15= 


Environment 




NOG^ 


group, 


field? 


group. 


Virgo 


group. 








off-center 




centre 


cluster 


near centre 



'NASA/IPAC Extragalactic Database 



"Lyon-Meudon Extragalactic Database 
'Third Reference Catalogue of Bright Galaxies 
iTruiilloetai] ( 120041) 
^This work 
iGiuricin eti5] ( l2000h 



rms scatters of points near the linear dependencies derived are 
about 0.2 A for all 4 indices under consid eration, within the ob- 
servational errors of IWorthev et al.l (Il994l) . 

To correct the index measurements for the stellar velocity 
dispersion which is usually substantial in the centres of early- 
type galaxies, we smoothed the spectrum of the K3III standard 
star HD 167042 by a set of Gaussians of various widths, and the 
derived dependencies of index corrections on cr were approxi- 
mated by 4h-order polynomials and applied to the measured in- 
dex values before their calibrations into the Lick system. Four 
galaxies of the five studied here have published ce ntral aperture 
inde x data obtained in the original Lick system (iTrager et al.l 
11998 ). In Fig. 2 we compare our fully calibrated central in- 
dex measurements aver aged between -2.2" and -1-2.2" for every 
galaxy with the data of ITrager et al.l (1 19981) taking in mind their 
aperture of 2" x 4". With 13 points of 16 within one sigma from 
the lines of equality, and with other 3 points within two sigma, 
the agreement seems to be good. 

To evaluate the internal accuracy of our index measurements, 
we made some simulations of the deepest galactic spectrum, 
the one for NGC 2634. The high S/N spectrum of the standard 
star HD 102328 was expanded along the slit with the normali- 
sation corresponding to the galaxy fluxes observed at different 
distances from the centre, the Gaussian smoothing was applied 
to imitate a radial velocity dispersion distribution, the sky level 
and the noise observed were added, and after that our standard 
reduction procedure was applied to this artificial spectrum. By 
measuring the indices, we assured that (i) difficulties related to 
the sky subtraction affect only the index measurements in the 
very external parts of the index profiles, beyond 2re. And (ii) 
the statistical accuracy of the indices estimated as the rms of the 
point-to-point scatter is 0.02-0.04 A over the range to 0.2 r^, 
0.03-0.07 A over the range 0.2 re to 0.5 re, 0.10-0.15 A over 
the range 0.5 to 1 re, and 0.21-0.27 A over the range 1 to 
2 re. As a visible illustration of the high internal accuracy of our 




SCORPIO SCORPIO 

Fig. 1. The correlation of the instrumental absorption-line in- 
dices obtained with the spectrograph SCORPIO and the stan- 
dard L ick indices for 22 stars from the list of IWorthev et al.l 
(11994 . The thin straight lines are bisectrices of the quadrants, 
and the dashed lines linear regressions fitted to the dependen- 
cies. The mean index error for the measurements of the stars by 
IWorthey et al.l (1 19941) are indicated for a single star in each plot; 
our errors are lower by one order 



index measurements. Fig. 3 compares two independent observa- 
tions of NGC 3348, made in 2005 and in 2006 with the different 
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SCORPIO SOORPIO 

Fig. 2. The comparison of the calibrated absorption-line indices 
obtained with the spectrograph SCORPIO for the central parts 
(1" X 3'.'5) of four galaxies and the s tandard Lick indic es for 
them from the aperture spectroscopy of lTrager et alj (119981) with 
the aperture of 2" x4". 
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slit orientations; the agreement of the indices is within 0.2 A up 
toR^ 30", or 1.36 re. 



4. The results 

4. 7 . Line-index profiies 

Figures 4 to 8 show the measured index profiles for all 5 galax- 
ies under consideration. Two halves of the profiles are symmetri- 
cally reflected with respect to the centre and plotted by different 
signs. The first thing that is evident from Figs. 4 to 8 is that the 
index radial variations cannot be fitted by a single linear law over 
the full radius range, in opposition to well-established tradition. 
With our high index accuracy, we clearly see breaks, or changes 
in slopes, of the observed dependencies. Below we give brief 
comments for every galaxy. 

NGC 2634, Fig. 4. This moderate-luminous elliptical galaxy 
was observed by us with the longest exposure of them all, so 
the observations are very deep. The index profiles are traced up 
to 3 Te, though we suspect that the sky subtraction uncertainty 
may affect their shapes at r > 2 r^. The emission lines are ab- 
solutely absent in the spectrum, and the H/3-index profile looks 
flat. The magnesium- and iron-index profiles decline with ra- 
dius. Their behaviour with respect to the logarithm of the radius 
may be characterised as a sequence of linear pieces with dif- 
ferent slopes. The points of slope breaks are close to 0.4 rg, or 
log(r/re) = -0.4, and, less surely, to 0.8 rg, or log(r/re) = -0.1. 
This break is seen in Mgb and Fe5335 profiles but not in the 
Fe 5270 profile. We take the value of for NGC 2634 from the 
work o fi Trujillo et al.l (l2004l) . where its photometric structure is 
described as a pure Sersic (non-core) one. The galaxy is a mem- 
ber of large group dominated by a giant Sb-galaxy NGC 2550, 
but it is located rather far from the centre of the group and has 



Fig. 3. Two independent observations of NGC 3348 in 2005 and 
in 2006 (small signs, filled and open) made with the orthog- 
onal slit or ie ntati ons. The aperture central measurements by 
iTrager et~arj (il998) are also shown by large filled circles with 
error bars. 

its own subgroup consisting of the dwarf NGC 2634A and of the 
late-type barred spiral NGC 2633. 

NGC 2810, Fig. 5. The observations are shallower reaching 
only 1.4 Tg. We determined the effective radius of the galaxy 
from our SCORPIO data in two ways: from the broad-band 
V-image that always precedes the spectral-mode exposure and 
from integration of the flux obtained in the spectral mode across 
the slit. The surfa ce brightness profiles are well-fitted by a de 
Vaucouleurs' law ( Ide Vaucouleurs|[T948l) . and both approaches 
have consistently given r^. - 14". The galaxy, though red and 
homogeneous photometrically, revealed strong emission lines in 
the spectrum, not only in the nucleus, but also over the full ra- 
dial extension. Due to large difference between the velocity dis- 
persion of stars and ionized gas in the nucleus, we succeeded 
in extracting a pure emission spectrum and estimated the ratio 
of the emission-line fluxes: Hj8/[Om] /i5007 = 0.85. Thi s ra- 
tio is within the value interval found bv lTrager et al.l (l2000al) for 
their sample of elliptical galaxies, 0.33-1.25, so we used it to 
correct the H/3 indices measured in the spectrum of NGC 2810 
for the emission through calculating EW([Om] A50Q1) over the 
full radial extension along the slit and multiplying it by 0.85. 
However, Fig. 5 presents the initial, uncorrected H/3 indices so 
the H/3 profile has a deep minimum in the centre of the galaxy. 
The Mgb index, very high in the centre, falls steeply down to 
log(r/re) ~ -0.1 and stays roughly constant beyond this ra- 
dius. The iron index decline is rather shallow over the full radius 
exte nsion. The el liptical galaxy NGC 28 10 is a ttributed neither 
bv iGarcial (Il993h nor by iGiuricin et al.l (l2000l) to any galactic 
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Fig. 4. The index profiles of NGC 2634. 
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Fig. 5. The index profiles of NGC 28 10. 



Fig. 6. The index profiles of NGC 3348. 



groups, so in our Table 1 we classify it formally as a field el- 
liptical. But according to HYPERLEDA, a lot of small late-type 
galaxies are seen not far from NGC 2810; some of them are even 
caught by the SCORPIO field of view. 

NGC 3348, Fig. 6. This galaxy is the most luminous in our 
sample, and it is the centre of its group. Truiillo et al. (200^ 
describe it as a core elliptical galaxy, and it is from this work 
that we take the effective radius value of 22". The spectrum 
demonstrates the weak emission line [Oiii] A5QQ7. As we were 
not able to extract the pure emission line HjS, we will correct the 
index H/? below for the emission in the statis t istical s ense, by 
following the recommendation of iTrager et al.l ( l2000al) : AHyS = 
0.6EW([Om] A5QQ1). Figure 6 presents an uncorrected H/3 pro- 
file, as well as magnesium- and iron-index profiles and com- 
bines the data for two different exposures of the galaxy, that of 
March 2005 and that of April 2006 made with the orthogonal 
orientations of the slit. The metal-line profiles fall smoothly to 
log(r/re) ~ -0.2. Beyond this radius , the Mgb and Fe5335 
profiles look flat, and the Fe5270 profile even rises. The final 
points of the profiles are around r « 1.3 rg. In April 2006, the 
faint galaxy at 2:5 to the east of NGC 3348, MCG -HI2-10-079, 
was also on the slit. It appeared to be the background early-type 
spiral galaxy, with the weak emission line [Om]/i5007 in the 
spectrum indicating ~ 10200 km s"' . Its nuclear Lick indices, 
H/3 = 1.38, Mgb = 3.32, Fe5270 = 3.03, and Fe5335 = 2.28, 
are typical of the bulge of the moderately luminous Sb galaxy. 

NGC 4339, Fig. 7. This early-type galaxy is a Virgo mem- 
ber It has a rather low luminosity for an elliptical and is clas- 
sified as either EO or SO by different authors. We analysed its 
surface brightness profile obtained by tracing the flux integrated 
across the slit and have seen that the surface brightness pro- 
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Fig. 7. The index profiles of NGC 4339. 



Fig. 8. The index profiles of IC 2179. 



file beyond r ^ 20" is well-fitted by an exponential, so be- 
yond r X 20" the stellar disk dominates. But within this radius 
we see a de Vaucouleurs' bulge with = 26", so this effec- 
tive radiu s is used to parametrise the index profiles in Fig. 7. 
ICaon. Ca paccioli. & D'Onofrio (1994) give an effective radius 
of 27"7 along the major and 25"8 along the minor axes, so in 
agreement with our measurements. Because of the bad weather, 
the exposure of NGC 4339 was short, and we only traced the 
index profiles for this galaxy to 0.5 rg. Nevertheless even these 
short profiles revealed the breaks of slopes, this time around 
log(r/re) ~ -0.7 — 0.6. The emission is absent in the galaxy, so 
the break in the H/3 profile at log(r/re) ^ -0.7 reflects the real 
change of the mean stellar age. 

IC 2179, Fig. 8. This elliptical galaxy, being the fainter 
neighbour of a giant spiral galaxy NGC 2347, has not been stud- 
ied in detail to now. In particular, the close bright star might 
preclude any attempts at surface photometry. We checked the 
brightness profile along the slit of the spectrograph and ensured 
that two halves of the profile diverge at r > 30". However, the 
surface brightness profile is very good inside this radius and can 
be fitted by a de Vaucouleurs' law (de Vaucouleurs 1948) with 
Te ~ 18"; if we take only the 'fainter' half of the profile, we 
obtain «; 15" in the radius range between 5" and 56" - this 
value is used to parametrise the profiles in Fig. 8. The metal-line 

profiles show two breaks - at log(r/re) ~ -0.3 0.2 and at 

log(r/re) ~ - -1-0.1; the H/3 index profile only reveals the latter 
break. The emission lines are absent in the spectrum, so the li/3 
index measurements are reliable. 



4.2. The magnesium-to-iron ratios 

Figure 9 presents the magnesium-iron diagrams for the in- 
dex profiles given in Figs. 4-8. We compare our observa- 
tional data with the evolutionary-sy nthesis models for SSP of 
iThomas. Maraston. & Benderl (|2003|) . which were calculated for 
several values of magnesium-to-iron abundance ratio, namely, 
the 0.5, 1, 2, and 3 solar ones ([Mg/Fe]=-0.3, 0.0, -hO.3, 
and -1-0.5). As it is well known elliptical galaxies have 
magnesium-overabundant stellar populations. Our results for 
the magnesium-to-iron ratio in the galaxies under consider- 
ation agree with all previous studies of the stellar popula- 
tion in elliptical galaxies. The mean values of [Mg/Fe] ra- 
tio change from -t-0.1 in NGC 4339 to +0.3 in NGC 3348, 
and our galaxies demonstrate clear dependence of the mag- 
nesium overabundanc e on the st e llar ve locity dispersion. 
As already noted by iTrager et al.l ( 12000 ^. the more mas- 
sive elliptical galaxies have higher magnesium-to-iron ra- 
tio. Within individual galaxies, the [Mg/Fe] ratios remain 
rou ghly constant along t h e rad i us as noted earlier by 
Wo rthev. Faber. & Gonzalez! d 19921), iDavies. Sadler. & Peletiej 
(11993*). and 'F isher. Franx. & Ilhngworthl (119951). The exception 
is NGC 28 10 where the centre has almost [Mg/Fe] - +0.5 - 
too high for its stellar velocity dispersion - and toward the outer 
part this ratio falls to its normal value, about +0.3. This elliptical 
galaxy has the most prominent emission lines in the spectrum; 
and as we see below, the young age of its central part is evidence 
of the secondary star formation burst perhaps provoked by a mi- 
nor merger. Less certain change of [Mg/Fe] from +0. 1 to +0.2 
can be noticed along the radius in NGC 4339. 
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Fig. 9. The magnesium-iron index diagrams for all the galaxies under consideration, with the Lick indices averaged over both sides 
of the slit taken alo ng the radius with the s teps as in Figs. 4-8 (circles connected by solid black lines, the nucleus marked as 'nuc'). 
The SSP models of lThomas. Maraston. & Bender (2003) for [Mg/Fe]=0.0, +0.3, and +0.5 (the last for NGC 2810 only) are plotted 
as a reference frame; the solid coloured lines represent stellar population models of equal ages; the metallicities for the models are 
+0.67, +0.35, 0.0, -0.33 from the top to bottom of each line. 
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4.3. The ages and metallicities along the radius 

In this subsection we present so-called SSP-equivalent ages and 
metallicities of the stellar populations determined by fitting the 
observed Lick indices H/3, Mgb, Fe5270, and Fe5335 by one- 
burst evolutionary synthesis models. The SSP models assume 
that all the stars within a stellar system have the same age and 
metallicity. In fact, the galaxies may be composite stellar sys- 
tems formed during several star formation bursts, with a differ- 
ent chemical composition in general. We therefore determined 
some a verage stellar ages an d metallicities. The recent simula- 
tions bv lSerra & Trageii (1200 6) have shown that the metallicities 
determined in this way are close to luminosity-weighted average 
metallicities, while the SSP-equivalent age is always biased to- 
ward the age of the youngest stars within the stellar system and 
is somewhere between the luminosity-weighted average age and 
that of the youngest stars. 

We consider together the ¥1/3 absorption-line index and the 
metal-Une combined one, 

[MgFe] = VMgb<Fe), (1) 

where 

Fe 5270 + Fe 5335 
(Fe) = , (2) 

to disentangle the age and metallicity effects and to deter- 
mine both characteristics simultaneously. The model calcula- 
tions show that, at the diagram H/3 vs [MgFe], the line families 
of equal metallicities and of equal ages (Figs. 10) have different 
slopes, so representing a somewhat curved coor dinate fra r ne. By 
using the model grid by Thomas, Maraston, & Bender! (l2003h 
for [Mg/Fe] - +0.3, we applied a 2D interpolation between the 
discrete model nodes to calculate the mean luminosity-weighted 
stellar ages and metallicities for each point of the profiles of 
Figs. 4-8. The H/3 indices in NGC 2810 and NGC 3348, where 
we see a significant emission contribution, are corrected for the 
emission by using the measurements of the equivalent width 
for the emission line [Oiii] /I5007: the additive correction is ob- 
tained by multiplying EW([0 Iii]/15007) by 0.85 for NGC 28 10 
(the individual coefficient estimated from the SCORPIO spec- 
tra) and by 0.6 for NGC 3348 (the statistical coefficient from 
iTrager et al.l(l2000ah ). 

To extract main tendencies of age and metallicity variations 
along the radius, we go on to average the indices, age, and metal- 
licity estimates over wider spatial bins near some key points of 
the profiles: re/8, re/5, re/2, r^, 1.3 re, and 2 rg. Here we calcu- 
late the errors of the mean indices and those of the mean ages 
and metallicities as simple point-to-point scatter divided by the 
number of measurements per bin. The results of this procedure 
are given in Table 3 and shown in Fig. 10 for each galaxy. 

Below we discuss the stellar metallicity and age variations 
along the radius in the individual galaxies. 

NGC 2634, Fig. 10a. This only galaxy in our sample has no 
stellar age gradient along the radius. The stellar population is 
old everywhere, 10-13 Gyr, and the mean metallicity changes 
from roughly -1-0.2 dex in the nucleus to less than -0.4 dex at 
r > rg. It may seem that the outer measurements of the in- 
dices are outside the parameter range of the models. But this 
discordance only re flects the need for more precision in tre ating 
low metallicities by iThomas. Maraston. & Bended (l2003l) (and 
by any other model sets as well). For comparison, we have plot- 
ted the index measurements for the Galactic globular clusters 
NGC 6624, NGC 6838, NG C 6356, and NGC 6539 with [Fe/H] 
from -0.35 to -0.66 from iBeaslev et al] (|2004|) : they are all near 



the ou ter measurements for NGC 2634 and outside the model 
grid of 'Thomas, Marast on. & Bended (12003 1). 

NGC 2810, Fig. 10b. The stellar population in the centre of 
this field elliptical is not very old, with a mean age of only 5-6 
Gyr. The metallicity of stars in the centre is very high, the high- 
est in our sample (though the galaxy is not the most massive 
and luminous one). Beyond the radius of 3" to 5", the stellar 
population is homogeneously old, about 10 Gyr, and the metal- 
licity is approximately solar. As we have also seen, there is a 
sharp increase in [Mg/Fe] ratio in the nucleus of NGC 28 10, so 
by taking a considerable amount of ionized gas in this galaxy 
into account, we may conclude that a minor merger has taken 
place some billion years ago to have provoked a short, intense 
nuclear star formation burst that increased the mean metallicity 
and decreased the mean age of stars. 

NGC 3348, Fig. 10c. The galaxy was classified by RSA as a 
regular EO, but it has now shown itself to be very unusual. The 
radial variations of the mean stellar age in this galaxy look like 
a mirror of NGC 2810: the nucleus is old, and the stars become 
younger and younger along the radius becoming 3^ Gyr old at 
r > rg. The metallicity gradient, though negative at r = 0-0.5 re, 
becomes positive in the radius range of 0.5 rg - 1.3 rg. Since 
NGC 3348 also shows extended emission lines in the spectrum, 
we may suggest that the galaxy has also experienced a minor 
merger But this time the satellite involved had a smaller amount 
of gas than that of NGC 28 10. As a result, no nuclear star forma- 
tion burst has been provoked, but the contribution of young and 
metal-rich stars of the merged galaxy is clearly seen in the outer 
part of NGC 3348. As evidence in favour of our hypothesis of a 
minor merger for NGC 3348, we demonstrate here the colour 
map of the galaxy (Fig. 11). We retrieved the HST/ACS im- 
ages of NGC 3348 obtained through two different filters, F435 W 
and F814W, on January 20, 2003, in the framework of the pro- 
gramme by W. E. Harris on globular cluster systems in early- 
type galaxies (Proposal ID 9427). By dividing one image by an- 
other and by taking 2.5 logarithms of the partial, we obtained 
the colour map, close to B - I, that is presented in Fig. 1 1 . One 
can see a lot of thin red filaments irregularly distributed over 
the whol e extension of the galax y. Evidently, this 'prototype 
of EO', as lSandage & Bedkd (Il99 4) characterized NGC 3348 in 
their Atlas, contains in fact a noticeable amount of diffuse dust 
and ionized gas coupled with this dust , which may be accreted 
rather recently. 

NGC 4339, Fig. lOd. Though less prominent than 
NGC 2810, NGC 4339 has also shown some signs of the 
circumnuclear young metal-rich stellar disk that may be a 
product of the secondary nuclear star formation burst. However, 
the gas is absent in this galaxy. 

IC2179, Fig. lOe. In IC2179 the nucleus is also younger 
than the outer part; but unless NGC 2810 and NGC 4339, 
IC2179 does not seem to have a young circumnuclear stellar 
disk producing a 'step' in age radial variations in the former 
galaxies. Instead we observe the gradual smooth age increase 
along the radius, from 8 Gyr to ~ 15 Gyr toward 2rg. The 
metallicity decrease along the radius is also very smooth, since 
IC2179 is the only galaxy in our sample where we do not see 
flattening of the metallicity gradient at r > 0.5 rg. 

5. Discussion 

The main point of our discussion concerns the metallicity gra- 
dients in elliptical galaxies, which are one of the signatures 
of their origin. As already mentioned, the monolithic dissi- 
pative collapse of a protogalactic cloud produces steep stel- 



M. Baes et al.: Metallicity and age gradients in round elliptical galaxies 

NGC 2634 NGC 2810 




2.5 



^ 2.0 



.5 - 



[Z/H] = -0.33 


















+0.67 
\ "^""^ 


















X T=5 

\^ T=8 


T=12 



[MgFe] 
NGC 3348 



.5 4 
[MgFe] 
NGG 4339 




IC 2179 



3.0 



2.5 



52, 

^ 2.0 



1 .5 



'1 

[Z/H] = -0.33 


'1 1 


1 ' 


\ \?*° 








+0.35 

Y-^^:L +0.67 
\ \ \ 









\^ T = S 

T=12 

T=llj 


2 


3 -1 
[K/qFe] 


5 



Fig. 10. The combined metal-line index vs diagrams for all the galaxies under consideration with the Lick indices averaged in 
wide bins near the key points along the radius - see the text and the Table 3 (circles with error bars, the nucleus being marked 
as 'nuc'). The SSP models of iThomas. Maraston. & Bender! (12003 ) for [Mg/Fe]= +0.3 are plotted as a reference frame; thin red 
and blue lines represent sequences of equal ages and metallicities , correspondingly. At the diagram for NGC 2634 several globular 
clusters of our Galaxy are also plotted as large red circles. 
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Table 3. The mean indices, ages, and metallicities along the radius in the galaxies of our sample 



galaxy 


log(r/re) 


Hj8 


Mgb 


<Fe) 


[Z/H] 


age [Gyr] 
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4.Dy + U.U/ 


Z.ao ± U.Uj 


+U.ZJ ± U.Uz 


11 o _L n A 
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4.j4 + U.Uj 


o An _i_ n n^ 
Z.W ± U.Uj 


+U. io ± U.U4 


1 o ^ _L n Q 
l/.j ± U.O 
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1.61 ±0.01 


3.89 + 0.06 


2.36 + 0.03 


-0.11 ± 0.02 


13.3 + 0.6 




0.0 


1.71+0.08 


3.25 + 0.08 


2.24 + 0.08 


-0.17 ±0.03 


9.5 ± 0.7 




0.3 


1.60 + 0.13 
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0.45 
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mean metallicity gradients obtained to date under the assump- 
tion of the constant stellar age along the radius from the colour 
and metal-line index gradients was about -0.2 dex per ra- 
dius dex, the opinion has formed that a major merger is typ- 
ical of a luminous elliptical galaxy to fla tten its me t allicity 
gradient. Refining this zero-order approach, 'Kobayashi' ('2004^, 
Forbes, Sanchez-Blazquez, & Proctor (2005), and Ogando et al] 
(t2005) all note a wide spread of metallicity gradients just for 
luminous (most massive) elliptical galaxies, thereby suggest- 
ing that there is a variety of evolutionary histories for elliptical 
galaxies, from quasi-monolithic collapse (or early gas-rich mul- 
tiple mergers) to 'dry' major mergers. 

In our work, we have disentangled the age and metallicity 
effects and measured both characteristics of the stellar popula- 
tions as they change along the radius. First of all, in most galax- 
ies the stellar age is not constant along the radius: its behaviour 
is quite various and does not depend on the galaxy's luminosity 
or environment. Secondly, when we have carefully analysed the 
metallicity profiles over the large radial base, to at least 1.3 r^, 
we find that the metallicity variations cannot be described by a 
single power law over the full radial extension in most galaxies. 
The break of the slope is evident at ~ 0.5 the outer parts show- 
ing the flatter metallicity gradients. We calculated the metallic- 
ity gradients for our 5 galaxies separately within r < 0.5 and 
in the radius range of 0.5 re - 1 r^, and the results are given in 
Table 4 with the galaxies ranked according to their stellar veloc- 
ity dispersion (mass). 




-Aa (arcsec) 

Fig. 11. The colour F435W-F814W (HST/ACS) map of 
NGC3348. 



lar m etallicity gradients along the r adius: from -0.5 (ICarlbergI 
[19841) to -1.0 dex per radius dex jKobavashil |2004|) . As the 
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Table 4. Metallicity gradients in our galaxies 
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galaxy 


d[Z/H]/dlogr 


d[Z/H]/dlogr 




(r < U.J r^) 


(0.5 Te < r < Te) 


NGC 3348 


-0.42 ± 0.05 


+0.20 ±0.17 


NGC 2810 


-0.88 ±0.17 


-0.47 ± 0.47 


IC2179 


-0.58 ±0.18 


-0.50 ± 0.40 


NGC 2634 


-0.52 ±0.16 


-0.20 ±0.17 


NGC 4339 


-0.82 ±0.52 





From Table 4 one can see that the metallicity gradients in the 
inner parts of the galaxies are rather steep, steeper than the -0.4 
dex per radius dex, and are consistent with the model of mono- 
lithic dissipative collapse. In the outer parts, the metallicity gra- 
dients are consistent with being zero; also taking the very inho- 
mogeneous mean age distributions at r > 0.5 re into account, we 
may suggest that the outer parts of the elliptical galaxies have 
been transformed in a secular manner through minor mergers, 
not by a major merger. This picture differs strongly from the 
conventional scenario of elliptical galaxy formation. 

An interesting question is the possible connection between 
metallicity and age variations. Earlier, Henrv & Worthev ( 199^ 
noted an anticorrelation between the metallicity and age gradi- 
ents found over the galaxy sample typically measured to 0.5 - 
1 rg. The slope of this relation appeared to be close to -|, which 
was a ratio of the sensitivities of the majority of meta l-line in- 
dices t o the age and meta l licity variations, AAge/AZ (IWorthevI 
11994 . iHenrv & WorthevI d 19991) propose that (i) all the differ- 
ences between the measured metallicity and age gradients in 
elliptical galaxies are due to random and systematic errors in 
index evaluation, producing the correlated scatter of the mea- 
sured metallicity and age gradients and (ii) that, in fact, all the 
elliptical galaxies have to have exactly d[Z/H]/dlogr - -0.25 
and dlogT/dlogr = -i-O.l. The accuracy of our measurements 
of the Lick indices along the radius of elliptical galaxies is 
higher than earlier results. We are therefore confident that the 
difference between the inner and outer metallicity gradients, say, 
in NGC 3348 or NGC 2634, is statistically significant. But in- 
deed, in our sample too, the metallicity and age gradients show 
a certain anticorrelation (if we exclude the very outer part of 
NGC2810), with a slope of -1.20 + 0.32. This is coincident 
with - 1 within the error Although this correlation may be real, 
its physical origin remains uncertain. A possible qualitative ex- 
planation is that additional star formation would simultaneously 
produce additional metals, while lowering the mean stellar age. 

Another important point concerning the observational signa- 
tures of the dominant mechanism for elliptical galaxy formation 
is correlation between the dynamical parameters, mass or escape 
velocity, and the mean abundance characteristics of the stellar 
population. Initially, the well-known correlation between lumi- 
nosity and the colou r of elliptica l galaxies was treated as mass- 
metallicity relation (lFabeii[l973b similar to what is known for 
other types of galaxies (e.g. lPilvu gin. Vflchez, & Contini 2004- 
But central-aperture spectral investigations of complete samples 
of nearby elliptical galaxies have given somewhat unexpected 
results: when disentangling the age-metallicity degeneracy and 
simultaneously measuring age, metallicity, and the [a/Fe] ratio, 
[Trageret al. (2000b) and recently Howell (2005) have noted a 
strong correlation between the [a/Fe] ratio and the stellar ve- 
locity dispersion, but the metallicity-cr, correlation appears to 
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Fig. 12. The correlation between the Mg b and (Fe) indices and 
metallicity and the local stellar velocity dispersion in three posi- 
tions along the radii in galaxies: in the nuclei, at r = 0.5 r^, and 
at r = r„. 



be weak. This fact provokes a revisi on of the formation models 
for elliptical galaxies. For example. IPipino & Matteuccil (|2004|) 
introduce an 'inverse wind' or a direct relation between the stel- 
lar velocity dispersion and star formation efficiency in every lo- 
cation within a proto-elliptical galaxy. High star-formation ef- 
ficiency in the densest regions means the shortest epoch of the 
main star formation, hence the highest [a/Fe] ratio that is con- 
sistent with the observational data; however, the physical base of 
this relation remains unclear. 

We have plotted our results concerning the relation between 
the metal-line indices and the mean stellar metallicity with the 
stellar velocity dispersion in Fig. 12; we have plotted not only 
the nuclei, but also the characteristics at r = 0.5 and at 
r = Tg. In the top plot, where we compare the Mgb index 
and the stellar velocity dis persion, we also indicate the mea n 
relations found earlier bv .Fisher. Franx. & IllingworthI (Il995h . 



relations round earl ier b y il^isher. j^ranx. & lllingwortnl 01^^31) . 
iTrager et alj (l2000bh . and lKuntschner et al. I (l2001h . The high ac- 
curacy of our data allows us to separate dependencies for the 
nuclei and for the off-nuclear regions. They go in parallel, with 
a significant vertical shift. Meanwhile, the relations published so 
far go much steeper than our observations and 'try' to involve 
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all the points that are in fact physically separated. The clearest 
physical picture is shown by the bottom plot, [Z/H] vs log cr. 
The points for r - 0.5 rg and r - are mixed; they represent 
good correlation, the rms scatter around the linear regression, 
[Z/H] - -1.54 + 0.69 logo-, is only 0.04 dex. In contrast, the 
points for the nuclei are certainly above this relation, and for 
the nuclei we do not see any correlation between the (nuclear) 
metallicity and the central stellar velocity dispersion. 

The separated position of the nuclear points in Fig. 12 forces 
us to suggest that the nuclear star formation in elliptical galaxies 
is never stopped by galactic wind, but instead proceeds quickly 
and efficiently toward the full exhaustion of the gas. The duration 
of star formation over the whole galactic body may be limited by 
the onset of a galactic wind due to the feedback from hot stars 
and supernovae and is governed by gravitation depending on the 
depth of the potential well. In any case, our results imply that 
all the physical correlations must be searched for only among 
the integrated characteristics of elliptical galaxies. Their stellar 
nuclei have their own evolution and are not representative of the 
galaxies as a whole. 



6. Conclusions 

By studying the Lick index gradients in 5 elliptical galaxies with 
moderate luminosity and located in different environments, we 
have carefully disentangled age-metallicity effects and measured 
radial variations of the mean luminosity-weighted stellar ages 
and metallicities up to 0.5 in one galaxy, up to 1.3 re in two 
galaxies and up to 2 rg in another two galaxies. We have found 
that the mean stellar age is constant along the radius only in one 
galaxy out of the 5. The other 4 galaxies demonstrate quite differ- 
ent mean stellar age behaviour: the outer parts are older than the 
centres in 3 cases and younger in one case. Among the former, 
we suspect the influence of the compact, rather young circumnu- 
clear disks in two galaxies. The metallicity gradients cannot be 
approximated by a single power law over the full extension of 
the radius in 4 galaxies out of 5. The inner metallicity gradients 
within 0.5 rg are all rather steep, steeper than -0.4 metallicity 
dex per radius dex, and are inconsistent with the origin of the 
elliptical galaxies by a major merger. The outer parts may rep- 
resent mostly the mix of stellar populations provided by minor 
mergers. 
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